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! Abstract 

We have done a series of two-dimensional hydrodynamic simulations of the rotational collapse of 

CO ' 

a supernova core in axisymmetry. We have employed a realistic equation of state (EOS) and taken 
into account electron captures and neutrino transport by the so-called leakage scheme. It is an 
important progress to apply the realistic EOS coupled with the microphysics to 2-D simulations for 
computing gravitational radiation in rotational core collapse. We have used the quadrupole formula 
to calculate the amplitudes and the waveforms of gravitational wave assuming the Newtonian 
gravity. From these computations, we have extended the conventional category of the gravitational 
waveforms. Our results have shown that the peak amplitudes of gravitational wave are mostly 
within the sensitivity range of the laser interferometers such as TAMA and first LIGO for a source 
at a distance of 10 kpc. Furthermore we have found that the amplitudes of the second peaks are 
within the detection limit of first LIGO for the source and first pointed out the importance of the 
detections, since they will give us the information as to the angular momentum distribution of 
evolved massive stars. 
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I. INTRODUCTION 



Asymmetric core collapse and supernova have been supposed to be one of the most plau- 
sible source of gravitational radiation for the long-baseline laser interferometers (GEO600, 
LIGO, TAMA, VIRGO) Q]. The detection of the gravitational signal is important not only 
for the direct confirmation of general relativity but also for the understanding of supernovae 
themselves because the gravitational wave is only a window that enables us to see directly 
the innermost part of an evolved star, where the angular momentum distribution and the 
equation of state are unknown. 

Observationally, the asymmetric aspects of the dynamics of supernovae are evident be- 
cause they are confirmed by many observations of SN 1987A and by the kick velocity 
of pulsars Q|. On the other hand, there is no consensus of the origin of asymmetry theoret- 
ically. However, provided the facts that the progenitors of collapse-driven supernovae are a 
rapid rotator on the main sequence [f| and that the recent theoretical studies suggest a fast 
rotating core prior to the collapse , rotation should play an important role as the origin of 
the asymmetric motions in the core collapse. It is noted that anisotropic neutrino radiations 
induced by rotation [8] could induce a jet-like explosion l| as suggested by the observations 
of SN 1987A. 

So far there have been works devoted to study the gravitational radiation in the rotational 
core collapse ^0 , 11, 12, 13, 14, 15, 1(| (see 01 for a review) . To say rigorously, reliable core 
collapse simulations should require the implementation of a realistic equation of state (EOS), 
an adequate treatment of microphysics (electron captures and other weak interactions) and 
neutrino transport, and a relativistic treatment of gravitation. However it is difficult to 
incorporate all of them at the same time. Therefore previous investigations have neglected 
or approximated the above requirements partially. So far most of the computations have 
oversimplified the microphysics. In addition, it has been reported that general relativity does 
not alter the significant features of gravitational radiation compared with those obtained in 
Newtonian approximation such as the range of gravitational wave amplitudes and frequencies 
1^ . This situation motivates us to employ a realistic EOS and treat microphysics adequately 
in the Newtonian gravity. In this paper, we will study the wave forms of gravitational 
radiation elaborately by performing the improved rotational core collapse simulations and 
will discuss what information can be extracted from the analysis. 



2 



We describe the numerical models in the next section. In the third section, we show the 
main numerical results. Conclusion is given in the last section. 



II. MODELS AND NUMERICAL METHODS 



Initial models 



We know little of angular momentum distributions in a core of evolved massive stars. 
Although it is supposed that some instabilities grow and transport angular momentum 
during the quasi-static evolutions, which mode prevails in what time scale is not understood 
very well at present. Therefore, we assume in this study the following two possible rotation 
laws. 

1. shellular rotation: 



R 
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Q{r) = Q ° x ^TM> (1) 

where Q(r) is an angular velocity, r is a radius, and Qq, Rq are model constants. 
2. cylindrical rotation: 

Q(X,Z) = Q x o X ° 2 9 ■ A Z ° (2) 

v ' ; x 2 + x 2 z* + z$' v ; 

where X and Z denote distances from the rotational axis and the equatorial plane and 
Xq, Z are model constants. The other parameters have the same meanings as above. 

Although recent theoretical studies (7[ give estimates for the angular velocity prior to 
the collapse, they are one-dimensional models with uncertainties and not the final answer. 
Hence we prefer a parametric approach in this paper. We have computed twelve models 
changing the combination of the total angular momentum, the rotation law, and the de- 
gree of differential rotation. The model parameters are presented in Table |l] The models 
are named after this combination, with the first letter , "S (Slow)", "M (Moderate)", "R 
(Rapid)", representing the initial T/|VK|init, the second letter, "S (Shellular), C (Cylindri- 
cal)", denoting the rotation law, and the third letter ,"L (Long), S(Short)", indicating the 
values of Rq and Xq, which represent the degree of differential rotation. The initial ratio 
of rotational energy to gravitational energy is designated as T l\W \\^\t- We have chosen 
0.25,0.5,1.5% for T/|W|i n it. The rotational progenitor model [18] corresponds to Model 
MSL (Moderate, Shellular rotation, Rq = 1000 km) in our simulations and we take this case 
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TABLE I: The Model Parameters. 



Model T/|W| init (%) Rotation Law R , X , Z xlO 8 (cm) Q (s^ 1 ) 



SSL 


0.25 


Shellular 


Ro = l 




2.8 


SSS 


0.25 


Shellular 


R = 0.1 




45.1 


SCL 


0.25 


Cylinder 


Xo = 1, Zo = 


1 


2.7 


SCS 


0.25 


Cylinder 


Xo = 0.1, Zo = 


= 1 


31.3 


MSL 


0.50 


Shellular 


Ro = l 




4.0 


MSS 


0.50 


Shellular 


Ro = 0.1 




63.4 


MCL 


0.50 


Cylinder 


Xo = 1, Zo = 


1 


3.8 


MCS 


0.50 


Cylinder 


Xo = 0.1, Zo = 


= 1 


44.4 


RSL 


1.50 


Shellular 


Ro = 1 




6.8 


RSS 


1.50 


Shellular 


Ro = 0.1 




112 


RCL 


1.50 


Cylinder 


Xo = 1, Zo = 


1 


6.6 


RCS 


1.50 


Cylinder 


Xo = 0.1, Zo = 


= 1 


76.8 



as the standard model. We have made precollapse models by taking a density and internal 
energy distribution from the spherically symmetric 15 M model by Woosley and Weaver 
(1995) and by adding the angular momentum according to the rotation laws stated above. 

B. Hydrodynamics 

The numerical method for hydrodynamic computations employed in this paper is based 
on the ZEUS-2D code |19| . The code is an Eulerian one based on the finite-difference 
method and employs an artificial viscosity of von Neumann and Richtmyer to capture shocks. 
The self-gravity is managed by solving the Poisson equation with the Incomplete Cholesky 
decomposition Conjugate Gradient (ICCG) method. Spherical coordinates are used and 
one quadrant of the meridian section is covered with 300 (r) x 50 (6) mesh points. The 
code is checked by standard tests such as the Sod shock-tube problem. We have made 
several major changes to the base code to include the microphysics. First we have added 
an equation for electron fraction to treat electron captures, which is solved separately. We 
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lave approximated electron captures and neutrino transport by the so-called leakage scheme 



20 



21 



22, 



. Second, we have in corp orated the tabulated equation of state (EOS) based 
on the relativistic mean field theory [24| instead of the ideal gas EOS assumed in the original 
code. It is noted that the implementation of the recent realistic EOS to 2-D simulations is 



an important progress beyond t 
the methods, see Kotake et al. 



le previous calculations. For a more detailed description of 



C. Gravitational wave signal 



We follow the methods by [10J, in order to compute the gravitational wave form. We 
will summarize it in the following for convenience. The amplitude of the gravitational wave 
h^ u = 9^u — Vfiv can be calculated by the quadrupole formula as follows: 

ft S T w = TrW^ ~ f »• < 3 > 

where i, j run from 1 to 3, t is the time, R is the distance from the source to the observer, 
the superscript "TT" means to take the transverse traceless part, and 1^ is the reduced 
quadrupole defined as 

Iij = J p{x){xiXj - ^x 2 5ij)d 3 x. (4) 

The second time derivative of the reduced quadrupole formula, which is difficult to be treated 
numerically, can be replaced by the spatial derivative by the equations of motion. In the 
case of axisymmetric collapse, the transverse traceless gravitational field is shown to have 
one independent component, hJ d T , and it is dependent solely on A^q . Then one derives for 
the component of h TT the following formula, 

where 9 is the polar angle and A^q is defined as 

^20 = ~T~7== / / drdfip[v r 2 (3fi 2 - 1) + v e 2 (2 - 3fi 2 ) - vJ - 6v r v g fiJl - fi 2 
c 4 V157T Jo Jo v 

-r<9 r $(3/i 2 - 1) + 3<%$ n^l-n*], (6) 

where d r = d/dr, dg = d/d8, p = cos 8 and $ is the gravitational potential. Since the 
gravitational wave is radiated most in the equatorial plane, the observer is assumed to be 
located in the plane in the following discussions. In addition, the source is assumed to be 
located at a distance of R = 10 kpc. 



TABLE II: Summary of important quantities for all models, is the time of bounce, /? max b is the 
maximum density at bounce, M{_ c b is the mass of inner core at bounce, T/|VF|fi na i is the final ratio 
of rotational energy to gravitational energy of the core, At is the duration time (FWHM) of the 
first burst, and |/i TT |max is the maximum amplitude of the first burst. Note that we speak of the 
inner core, where the matter falls subsonically, which corresponds to the unshocked region after 
core bounce. 



Model 


tb 


Pmaxb 


M;. cb 


T/\W\ &nal 


At 


\h TT \ 

| " | max 




(ms) 


(10 14 g cm" 3 ) 


(M e ) 


(%) 


(ms) 


(io- 20 ) 


SSL 


227.4 


2.95 


0.74 


4.87 


0.58 


1.03 


SSS 


227.3 


2.80 


0.75 


5.66 


0.54 


1.49 


SCL 


241.2 


2.87 


0.75 


5.09 


0.46 


0.93 


scs 


230.7 


2.84 


0.73 


6.46 


0.56 


2.00 


MSL 


242.7 


2.65 


0.75 


8.44 


0.72 


1.58 


MSS 


242.2 


2.15 


0.87 


9.05 


0.57 


2.03 


MCL 


241.9 


2.85 


0.76 


8.38 


0.70 


1.53 


MCS 


245.4 


1.45 


0.91 


9.94 


0.51 


2.87 


RSL 


338.7 


1.21 


0.94 


14.6 


2.89 


0.48 


RSS 


328.6 


0.53 


0.96 


13.3 


2.80 


0.78 


RCL 


338.4 


1.36 


0.95 


14.5 


4.47 


0.42 


RCS 


326.7 


0.18 


1.11 


12.3 


1.91 


1.41 



III. NUMERICAL RESULTS 

1. The Properties of the Waveform 

We first show the general properties of the waveform with collapse dynamics. We choose 
model MSL (standard) as a representative model. For later convenience, the values of several 
important quantities are summarized in Table |TI] The time evolution of the amplitude of 
gravitational wave and the central density near core bounce are shown in Figure As 
the inner core shrinks, the central density increases and a core bounce occurs when the 
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central density reaches its peak at % = 243 msec with 2.65 10 14 g cm -3 . At this time, the 
absolute value of the amplitude becomes maximum. After the core bounce, the core slightly 
re-expands and oscillates around its equilibrium. As a result, the gravitational wave shows 
several small bursts and begins to decay. These gross properties are common to all other 
models. However, there exist some important differences when we compare them in more 
detail. We will discuss the differences in the following. 

In Figure El the time evolution of the amplitude of gravitational wave and the central 
density near core bounce for model MCS are given. For the model, it is noted that the initial 
rotation law is cylindrical with strong differential rotation. By comparing the left panel of 
Figure Q to El the oscillation period of the inner core for model MCS is clearly longer than 
for model MSL. In other words, the pronounced peaks can be seen distinctively in this case. 
This is because the central density becomes more smaller after the distinct bursts (compare 
the right panels of Figure ^ to EJ . This effect increases with the initial angular momentum 
(compare the left panel of Figure El t° El)- It is also found that the signs of the values of 
the second peaks are negative for model MCS, on the other hand, positive for model MSL 
(compare the left panels of Figure El to QJ. Note that we will speak of the second peak 
where the absolute amplitude is second largest. The above characteristics are common to 
the models for strongly differential rotation with cylindrical rotation law (see Table ITTT^ . A 
specific feature in the waveform is found for models RSL and RSS in which there exist the 
small peak or shoulder before the peak burst (see the right panel of Figure Ell- 
Next we will compare our results of the waveforms with those by Zwerger et al. 12] who 
categorized the shapes of waveforms into three district classes. They used a polytropic EOS 
to express the pressure from the degenerate leptons as P oc p 1 . They reduced the adiabatic 
index, 7, from 1.325 to 1.28 below the nuclear density regime in order to approximate the 
related microphysics and neutrino transport. They assumed the cylindrical rotation law 
for all the models and varied both the degree of differential rotation and the initial angular 
momentum. Our models except for the strongly differential rotation with cylindrical rotation 
law correspond to so-called the type I in their nomenclature. On the other hand, our models 
for strongly differential rotation with cylindrical (not shellular) rotation law correspond to 
type II. It should be noted that type II signals were limited to occur only for rapid, strongly 
differential rotation with cylindrical rotation law in their models. Furthermore in this study, 
we find that type II does not occur for models with shellular rotation law regardless of the 
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TABLE III: Some characteristic quantities for the waveform analysis. The names of the models 
whose initial rotation law is cylindrical with strong differential rotation are written in bold letters. 
Tq SC and T**, the first and second oscillation period of the inner core, respectively. h^ ond is the 
amplitude of gravitational wave at the second peak. 



Model Ti c 
[ms] 


osc 

[ms] 


second 

[io- 20 ] 


SSL 


1.6 


2.4 


0.82 


SSS 


1.5 


2.2 


0.94 


SCL 


1.6 


2.3 


0.75 


SCS 


1.8 


2.6 


-0.57 


MSL 


1.8 


1.2 


1.23 


MSS 


2.2 


2.4 


0.98 


MCL 


1.7 


1.1 


1.20 


MCS 


3.1 


2.6 


-0.79 


RSL 


1.0 


3.1 


0.32 


RSS 


1.0 


7.8 


0.20 


RCL 


9.0 


2.9 


0.25 


RCS 


10.7 


8.6 


-0.49 



degree of differential rotation and the initial rotation rate, on the other hand, does occur 
regardless of the initial rotation rate in case of strong differential rotation with cylindrical 
rotation law. Finally type III occurred only when the core collapses very rapidly (7 = 1.28) 
in their calculations. There are no models which correspond to the type in our calculations. 
This is because our realistic EOS is not so soft in the corresponding density regime. 

2. Maximum Amplitude and Second Peak 

We will first discuss the relation between the maximum amplitudes of gravitational wave 
and the initial T/IW 7 !. The maximum values and the related quantities are given in Table 
H] and |n] From Figure EJ it is found that the amplitude is largest for moderate initial 
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rotation rate (i.e., T/|VK|hut — 0.5%) when one fixes an initial rotation law and a degree 
of differential rotation. This is understood as follows. The amplitude of gravitational wave 
is roughly proportional to the inverse square of the typical dynamical scale, t dyn (e.g., Eq. 
(J21) )• Since tdyn is proportional to the inverse root of the central density p, the amplitude 
is proportional to the density. As a result, the amplitude becomes smaller as the initial 
rotation rates become larger because the density decreases for large initial rotation rate (see 
Table HTl and Figure EJ). On the other hand, the amplitude is proportional to the value of the 
quadrupole moment, which becomes large in turn as the total angular momentum increases. 
This is because the stronger centrifugal forces make not only the mass of the inner core 
larger (see Table ITTj) . but also deforms the inner core. The amplitude is determined by the 
competition of these factors. As a result, the amplitude becomes extrema for moderate 
initial rotation rates in our calculation. 

Next we will discuss the values of maximum amplitude between our results and those by 
other groups. The values of maximum amplitude for all our models range from 5 x 10 -21 < 
h TT < 3 x 10~ 20 , which is almost the same as the results by Zwerger et al. [12] and 
Monchmeyer et al. [10]. On the other hand, the values of the standard models by Yamada 

II] are about an order of magnitude lower. This is understood as follows. Yamada 

III ] employed a parametric EOS by which the effects of microphysical and transport 



et al. 
et al. 



processes were assumed to be expressed. Since they found that the maximum amplitude is 
most sensitive to the adiabatic index at the subnuclear density by their parametric surveys, 
we pay attention to this. The comparison of the effective adiabatic index between our EOS 
and their EOS is shown in Figure El Compared to their EOS, our realistic EOS is rather 
soft for the subnuclear density regime. At the regime, our realistic EOS can express the 
softening of EOS by the effect that the nuclear interactions become attractive. Due to this 
effect, the inner core can shrink more compact at core bounce, which results in the larger 
maximum amplitude. It is naturally suggested that we may get the information about the 
subnuclear matter if we can detect the gravitational wave from the rotationally collapsing 
cores. We hope it can be realized in the near future since the maximum amplitudes for our 
models are mostly within the detection limit for TAMA and first LIGO which are now in 
operation if a source is located at a distance of 10 kpc (see Figured). 

As pointed earlier, we find that the signs of the values of the second peaks are negative 
for the models for strongly differential rotation with cylindrical rotation law and positive for 
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the others (see Table ITTTj) . The absolute amplitudes of the second peak are also presented 
in Figure [7| As shown, they are within the detection limit of first LIGO for a source at 
a distance of 10 kpc. In addition, it seems quite possible for the detectors in the next 
generation such as advanced LIGO and LCGT to detect the difference. Therefore if we can 
find the difference of the signs of the second peaks by observations of gravitational wave, we 
will obtain the information about the angular momentum distribution of evolved massive 
stars. Since there is no way except for the detection of the gravitational wave to obtain the 
information, it is of great importance to detect the second peaks in the future. 



3. Secular instability 

In some models (Models RSL and RCS) with large initial rotation rate, the final rotation 
rate exceeded the critical value (see Table HT]) . where MacLaurin spheroids become secularly 
unstable against tri-axial perturbations (T/|VK| > T/|VK| seq = 13.75%). Rampp et al. [L3j 
reported that no considerable enhancement of gravitational radiation due to the growth of 
secular instabilities was found within a time scale of several 10 ms after core bounce. Since 
our calculations are within the time scale, the axial symmetry assumed in this work may be 
justified. 



IV. CONCLUSION 



We have done a series of two-dimensional hydrodynamic simulations of the rotational 
collapse of a supernova core and calculated gravitational waveforms using the quadrupole 
formula. We have employed a realistic EOS and taken into account electron captures and 
neutrino transport in an approximated method. We have found the following: 

1. The peak amplitudes of gravitational wave obtained in this study are mostly within 
the detection limits of the detectors of TAMA and first LIGO which are now in operation if 
a source is located at a distance of 10 kpc. In addition, the peak amplitude becomes extrema 
for the models whose initial rotation rate is moderate. 

2. The waveforms are categorized into the criteria by Zwerger et al. Q|. In addition, 
we further find that type II does not occur for models with shellular rotation law, on the 
other hand, does occur regardless of the initial rotation rate in case of strong differential 
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rotation with cylindrical rotation law, and that the type III does not occur if a realistic EOS 
is employed. 

3. At the subnuclear density regime, our realistic EOS can express the softening of EOS 
by the effect that the nuclear interactions become attractive. Therefore the inner core can 
shrink more compact at core bounce than the other work in which EOS is expressed 
in a parametric manner. Subsequently, this results in the larger maximum amplitude. It 
follows that we may get the information about the subnuclear matter if we can detect the 
gravitational wave from the rotationally collapsing cores. 

4. The signs of the values of the second peaks are negative for the models with strong 
differential rotation with the cylindrical rotation law, on the contrary, positive for the other 
models. The absolute amplitudes of the second peaks are within the detection limit of first 
LIGO. Therefore if we can detect the signs of the second peaks, it will give us the information 
as to the angular momentum distribution of massive evolved stars when a supernova occurs 
at our galactic center. 

As stated earlier, the detection of gravitational wave is likely for the models whose initial 
rotation rate is moderate. According to the study of rotational core collapse by Kotake et 
al. the anisotropic neutrino radiation is well induced by such a rotation rate. Noting 
that the anisotropic neutrino radiation can induce the globally asymmetric explosion jj 
the detection of gravitational wave will be a good tool to help understand the explosion 
mechanism itself. 
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240 245 250 255 260 265 240 242 244 246 248 250 252 

Time [ms] Time [ms] 

FIG. 1: Time evolution of the amplitude of gravitational wave (left panel) and central density near 
core bounce (right panel) for model MSL. 




240 245 250 255 260 265 240 242 244 246 248 250 252 254 

Time [ms] Time [ms] 

FIG. 2: Same as Figure □ but for model MCS. 
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FIG. 3: Waveforms for models RCS (left panel) and RSL (right panel). For the left panel, very 
district peaks are shown. For the right panel, a small peak and shoulder after the first burst can 
be seen. 




T/|W| init [%] 



FIG. 4: Relation between T/|W|init an d the peak amplitude |/i TT | m ax for all the models. "CS, SS, 
SL, CL seq" in the figure represent the model sequences whose names are taken from the second 
and third letters in Tabled respectively. Note that the second and third letters mean the rotation 
law and the degree of differential rotation, respectively. 
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3.5 



CS seq 
SS seq 
SL seq 
CL seq 




T/|W| init [%] 



FIG. 5: Relation between T/|W|i n it and the central density at core bounce (10 14 g cm -3 ) for all 
the models. The meaning of the label of lines is the same as Figure 01 It is found that the central 
density at the core bounce becomes smaller as the initial angular momentum becomes large. 




1.2 ' 1 1 1 1 1 1 

13 13.2 13.4 13.6 13.8 14 14.2 

log(Density [g/cm 3 ]) 



FIG. 6: Relation between the density and the effective adiabatic index 7 at the subnuclear density. 
MSL represents the relation for our standard model. On the other hand, the index employed in 
the study of Yamada et al. is represented as YS. 
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FIG. 7: Detection limits of TAMA [2!|, first LIGO |2f|, advanced LIGO |23], and LCGT 
the expected amplitudes from numerical simulations. The open squares represent the maximum 
amplitudes for all the models. On the other hand, the pluses and the closed squares represent the 
amplitudes of second peaks for models with strong differential rotation with cylindrical rotation law 
and for the other models, respectively. We estimate the characteristic frequencies by the inverse 
of duration periods of the corresponding peaks. Note that the source is assumed to be located at 
the distance of 10 kpc. 
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